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Abstract: This study was to characterize dynamic source strength changes estimated from high-density scalp 
electroencephalogram (EEG) at different phases of a submaximal voluntary muscle contraction. Eight healthy volunteers 
performed isometric handgrip contractions of the right arm at 20% maximal intensity. Signals of the handgrip force, 
electromyography (EMG) from the finger flexor and extensor muscles and 64-channel EEG were acquired 
simultaneously. Sources of the EEG were analyzed at 19 time points across preparation, execution and sustaining phases 
of the handgrip. A 3-layer boundary element model (BEM) based on the MNI (Montreal Neurological Institute) brain 
MRI was used to overlay the sources. A distributed current density model, LORETA LI norm method was applied to the 
data that had been processed by independent component analysis (ICA). Statistical analysis based on a mixed-effects 
polynomial regression model showed a significant and consistent time-dependent non-linear source strength change 
pattern in different phases of the handgrip. The source strength increased at the preparation phase, peaked at the force 
onset time and decreased in the sustaining phase. There was no significant difference in the changing pattern of the source 
strength among Brodmann's areas 1, 2, 3, 4, and 6. These results show, for the first time, a high time resolution 
increasing-and-decreasing pattern of activation among the sensorimotor regions with the highest activity occurs at the 
muscle activity onset. The similarity in the source strength time courses among the cortical centers examined suggests a 
synchronized parallel function in controlling the motor activity. 

Keywords: EEG source, Current density reconstruction, Electroencephalography (EEG), Brain, Handgrip force, Voluntary 
muscle contraction. 



INTRODUCTION 

The time course of brain activation during voluntary 
muscle contraction was first reported by the studies invol- 
ving single-cell recordings in the brain of non-human pri- 
mates [1-3]. These studies reported varying activation levels 
of motor cortical neurons among various phases (e.g., pre- 
paration vs. execution) during a motor control process. 
Although the observations provide valuable information 
describing single cell activities with time in animals in 
controlling a voluntary motor action, very little is known 
about activation time course of individual brain regions in 
humans in millisecond resolution. 

Time courses of human cortical field activation are typi- 
cally studied by neuroimaging involving relatively long- 
duration tasks. For example, alterations of activation of 
various brain areas were investigated during muscle fatigue 
lasting for many seconds using functional magnetic reso- 
nance imaging (fMRI) [4, 5] and positron emission tomo- 
graphy (PET) [6]. These studies demonstrate time- or 
fatigue-dependent alterations of cortical activities across a 
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relatively long time but cannot identify potential shifts of the 
activity level from one stage (e.g., planning) to another (e.g., 
execution) in a motor control process given the relatively 
poor temporal resolution of current-available neural imaging 
techniques. Although scalp electroencephalography (EEG) 
[7] or magnetoencephalography (MEG) [8] signals have 
millisecond resolution and are frequently recorded to 
describe control strategies of the brain, the source(s) of the 
signals or the involving cortical fields are usually unknown 
because of relatively poor spatial resolution of these 
methods. To improve spatial resolution and maintain high 
temporal resolution of scalp EEG or MEG data, several 
inverse solution methods to estimate the signal sources have 
been developed [9]. These inverse methods aim to de-blur 
the head volume conductor effect on the signal recorded 
from the scalp to identify the true electrophysiological 
sources in the cortex. With detailed temporal and spatial 
information, the inverse method could serve as a promising 
technique to explore activation time course of a cortical 
region or regions in movement control. 

Commonly used inverse models include equivalent- 
current-dipole, distributed current density, beamforming, etc. 
Equivalent-current-dipole model relies on an ad hoc assump- 
tion of a reasonable number of dipoles, their known locations 
and/or orientations. Distributed current density model does 
not need prior knowledge of specific source number and 
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locations. Among the several available current-density 
estimation techniques, the low-resolution electromagnetic 
tomography (LORETA) has been shown to have better 
source localization ability than other methods [10, 11]. The 
LORETA technique has been used in investigating abnor- 
malities of P300 in schizophrenia [12, 13], analyzing wide- 
spread epileptiform activities in epilepsy [14, 15], and exp- 
loring cognitive function of the brain in healthy individuals 
[16]. 

The purpose of this study was to use LORETA source- 
localization technique to estimate sources of scalp EEG 
signals (64 channels) recorded before and during a handgrip 
motor task and determine the time-dependent source strength 
(activity level) changes across preparation, execution and 
sustaining phases of the task in different cortical sensori- 
motor areas. It was hypothesized that neural activity in the 
cortical regions would vary depending on the control phases 
of the handgrip task and the changes could be detected by 
the current density reconstruction technique. 

MATERIALS AND METHODOLOGY 

Subjects and Motor Task 

Eight right-handed subjects (4 men and 4 woman, age = 
31 ± 2 yrs) participated in the study, none had neurological 
or other disorders at the time of study. The experimental 
procedures were approved by the Institutional Review Board 
at the Cleveland Clinic. All subjects gave informed consent 
prior to the participation. 

Subjects were seated comfortably in an experimental 
chair in an electrically shielded data-recording room. During 
the experiment, subjects performed 20 isometric handgrip 
contractions of the right arm at 20% maximal voluntary 
contraction (MVC) level. An oscilloscope was placed in 
front of the subjects and the target was displayed on the 
screen. In each trial, they steadily increased the handgrip 
force to the target (20% MVC) level and matched the exerted 
force as closely as possible with the target for about 3 s. 

Special precautions were taken to minimize noises in the 
data during the experiments. The subjects were required not 
to move the head and body, blink the eyes, bite the teeth, 
contract the facial/neck muscles, or tense the body while 
performing the low-intensity handgrips. Muscle fatigue was 
minimal by allowing sufficient rest between adjacent trials. 
The length of the inter-trial rest period was controlled by the 
subjects themselves, i.e., they started the next trial when they 
felt comfortable and completely rested. The inter-trial rest 
period was, on average about 20s. 

Data Recording 

Force 

The handgrip force was recorded by a pressure 
transducer (EPX-N1 250 PSI, Entran Devices, Inc., Fairfield, 
NJ) interfaced with distill water in a bottle-like device 
through a nylon tube [17]. The output of the transducer was 
directed to the amplifier and then to an input channel of a 
Micro 1401 data acquisition system (version 3.05, 
Cambridge Electronic Design, Ltd., Cambridge, UK), which 
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transferred the voltage data to a laptop computer. The 
sampling rate for digitizing the force data was 100 Hz. 

EEG 

A 64-channel NeuroScan EEG system (version 4.2, 
NeuroScan, El Paso, Texas, USA) was used to acquire EEG 
signals from the scalp. A Quik-Cap elastic nylon cap that 
held 64 electrodes was placed on the scalp for EEG data 
recording. Conducting gel was injected into each electrode to 
connect the recording surface of the electrode with the scalp. 
Impedance was controlled below 10,000 ohms. One elec- 
trode was used to record the force, which was used for the 
purpose of synchronizing the EEG with muscle activities. 
All the EEG electrodes were referenced to the common 
linked earlobes and the signals were amplified (X75,000), 
low-pass filtered (0.05-50 Hz), and digitized (2,000 
sample/s). 

Data Analysis 

The force data were processed using house-coded 
programs within the Spike2 software package (Cambridge 
Electronic Design, Ltd., Cambridge, UK). The force onset 
time was set at the 10% handgrip MVC force rising from the 
baseline in each trial, which was used for the purpose of 
synchronizing the force, and EEG signals. The actual force 
was measured by averaging the data points in each trial when 
the force was steady (rising and falling portions were 
excluded) and then averaging across the trials before 
performing the group averaging. 

The raw EEG data were visually inspected and trials with 
artifacts due to eye blinks or head movements were excluded 
(on average, 2+2 trials were removed in each subject). Both 
EEG data preprocessing and LORETA current density 
estimation were performed using Curry software package 
(version 4.5, Neuro Scan Labs, Virginia, USA). 

Estimation of Current Densities Using LORETA 

For each subject the EEG signals were aligned with the 
force onset and then averaged across trials. The averaged 
EEG data were baseline corrected using baseline data from - 
3000ms to -2500ms (time 0 indicates force onset). 
Subsequently, an independent component analysis (ICA) 
was applied to the data. Only the main components (signal- 
to-noise ratio > 1) were chosen for the source reconstruction. 
A 3-layer (conductivities of the scalp and brain: 0.033S/m, 
and the skull: 0.0042S/m) triangle-node boundary element 
model (BEM) [18, 19] based on the MNI (Montreal 
Neurological Institute) brain MRI was used to overlay the 
sources. Distributed current density model (LORETA) with 
LI norm method was applied to the ICA-preprocessed data. 
In addition, the sources were constrained to the reconstructed 
layer of the folded cortex with 6926 nodes using a rotating 
model [Curry user guide, 1999, 11]. 

Current Density Data Analysis - Post Processing 

Source localization was analyzed at 21 different time 
points (-2200ms, -2100ms, -2000ms,- 1000ms, -900ms, 
-800ms, -700ms, -600ms, -500ms, -400ms, -300ms, -200ms, 
-100ms, 0ms, 100ms, 200ms, 500ms, 800ms, 1000ms, 
2000ms, 2500ms and 3000ms) throughout the planning, 
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Fig. (1). Reconstructed normalized and averaged cortical currents at 20% maximal voluntary contraction force of 8 subjects at eight time 
points from early preparation to holding phases of the handgrip task. The size of the red dots is proportional to the local current density. Only 
currents in motor related cortical area (left Ml, left SI, PM&SMA, right Ml, right SI) are shown here. The current density exhibited non- 
linear alterations throughout the motor control process. 



execution, and sustaining phases of the handgrip muscle con- 
traction in each subject. The current density at each location 
was averaged across subjects. Since each source needed to 
be identified at its anatomical location under Talairach 
coordinate for further analysis, several steps were taken to 
transform Curry V4. 5 coordinate to Talairach coordinate. 
First, transformation of the SPM99/ MNI (X,Y,Z) coordi- 
nates from the Curry coordinates (x,y,z) was obtained as 
follows (the MNI image dataset has a 1.8-mm voxel size in 
Curry V4.5, while the MNI brain originally had a voxel size 
of 2 mm): 

X = (120mm-x)*2/1.8 

Y = (102mm-y)*2/1.8 

Z = (z- 100 mm) * 2 / 1.8 

Second, transformation of the SPM99/ MNI coordinates 
(X, Y, Z) to the Talairach coordinates (X 1 , Y', Z') was per- 
formed: (http://www.mrc-cbu.cam.ac.uk/Imaging/Common/ 
mnispace.shtml) 

Above the anterior commisure (AC) (Z >= 0): 

X'= 0.9900X 
Y'= 0.9688Y +0.0460Z 
Z'= -0.0485Y +0.9189Z 
Below the AC (Z < 0): 

X'= 0.9900X 

Y'= 0.9688Y +0.0420Z 

Z'= -0.0485Y +0.8390Z 

Third, after establishing the position of each source on 
the Talairach coordinate, the anatomical label was obtained 



through the Talairach Daemon search for each source 
location. Thus, all the current densities in the Brodmann's 
area 6 (premotor [PM] and supplementary motor area 
[SMA]), area 4 (primary motor cortex [Ml]) and areas 1, 2 
and 3 (primary sensory cortex [SI]) were identified and 
vector-averaged. (We chose the three Brodmann's areas 
based on their well known role in regulating muscle activi- 
ties.) The overall averaged current densities of these three 
major sensorimotor areas at time = -2200ms, -2100ms and 
-2000 ms were used as the normalization factor for each 
subject. The purpose of this normalization was to reduce the 
confounding effects of possible wide variation of global 
neuronal activity among subjects. Finally, the current density 
at each time point was normalized by this normalization 
factor for further statistical analysis. 

Statistical Analysis 

This study represented a randomized nested design, in 
which each observed curve was classified according to the 
factor "subject" and factor "area" of the subject. The nature 
of the experimental data was such that the curves of current 
density (CD) by various Brodmann's areas can be different, 
but the difference is not expected to be systematic in terms 
of area. Thus "subject" and "area" are considered to be nes- 
ted classification factors. Since eight subjects were randomly 
chosen, random- effects terms needed to be associated with 
the "subject" factor and with the "area" factor nested within 
a subject. 

To model the nonlinear pattern of the source strength 
variation, a mixed-effects polynomial regression model [20] 
for the CD data with respect to the "time" covariate and 
"area" factor was considered. We used random-effects for 
the intercept, linear term and quadratic term at the subject 
level and a single random effect for the intercept at the area 



Cortical Activation in Handgrip Contraction 



The Open Neuroimaging Journal, 2011, Volume 5 235 



within subject level. This setting allowed the overall pattern 
to vary between subjects not only in the CD peak location 
but also the curvature at the peak. The difference between 
areas for the same subject is a shift in the intercept. Akaike 
information criterion (AIC) [21] was used for modeling 
selection and to determine the order of the polynomial 
function in the regression model. It turned out that the fourth 
order polynomial function was adequate to model the 
nonlinear pattern of the source strength over time. 

If CD of the jth area was written on the ith subject at the 
kth time point as y^, i=l,..., 8; j=l,..., 5; k=l, K, the 
final model being fit can be expressed as 

y.jt = Po + Pih + PA + PA + PA + area j + b l0 + b u t k + b i2 t 2 t + b tj + e ijt , 

where bj,o, bjj, b;,2 are random effects terms (following 
normal distributions with unknown variance parameters) 
associated with subject factor, by is the random effects term 
associated with area factor nested within subject factor. Ep 
are random noise in the model. Maximum likelihood (ML) 
approach with a ridge-stabilized Newton-Raphson algorithm 
approach [22] was used to fit our specified model. ANOVA 
F tests [22] were constructed based on the mixed models to 
test "area" effect and time trends at the 0.05 significance 
level. 



RESULTS 

Fig. (1) shows the normalized and averaged (based on the 
8 subjects) current density (CD) reconstruction maps in the 
primary motor cortex (Ml), primary sensory cortex (SI), 
premotor cortex (PM), and supplementary motor area (SMA) 
at different time points before and after the force onset. It is 
clear that the EEG source strength (CD) in the cortical areas 
progressively increased from the early preparation phase to 
the movement execution phase, especially for the Ml and SI 
regions. The CD in the SMA increased earlier during the 
preparation and localized predominantly in the left frontal 
lobe. 

Fig. (2) shows the normalized source strength calculated 
at all time points in the left Ml, left SI, the PM and SMA 
(calculated bilaterally), right Ml, and right SI . The CD in all 
the areas began to rise early from ~ 1000ms before the force 
onset, peaked at the force onset, and then drifted downward 
afterwards. During the entire time course, the left Ml and SI 
exhibited greater source strength compared with the other 
areas and source strength of the Ml and SI showed a 
tendency to rise again 500ms after force onset. 

The mixed model provides us with the flexibility of 
modeling not only the mean structures of our data but their 
variances and covariances as well. Specified random effects 
at different levels allowed us to fit a specific curve for each 
subject at each area within a unified model. Here we used 
normalized time when the model was fitted (i.e. 0<t<l). 
Analysis of covariance had been applied to the mixed model 
for "source strength" as it relates to both a classification 
factor (area) and to a continuous covariate (time). To make it 
easier to assess the differences between areas, we used "left 
Ml" representing an "overall mean" or reference level and 
other four areas representing changes to the reference level. 
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Fig. (2). Averaged and normalized current density plotted as a 
function of time in the motor related cortical area (left Ml, left SI, 
PM&SMA, right Ml, right SI) in 8 the subjects at 20% maximal 
voluntary handgrip force. Time '0' is the force onset. Negative time 
values indicate time points before force onset. 

Table 1 displays the statistical results of analysis of 
covariance based on the mixed-effects polynomial regression 
model. The "Time" variable was normalized into the interval 
[0,1] in the model fitting. The p-values were presented for 
global fixed effects over all subjects. It is noted that all terms 
of polynomial coefficients were significant. Since the factor 
"area" was a categorical variable that had five categories, the 
"left motor" area was set as the baseline and its estimate and 
standard error were zeros. The ANOVA test for the factor of 
"area" was not significant with a p-value = 0.7255. These 
results indicated that there was no statistical difference in 
source strength among the areas but the strength changed 
significantly over time in a nonlinear function (Fig. 3). Fig. 
(3) shows the original data and fitted curves for each area 
within each subject from the mixed model. The model 
included random effects from both subject level and area 
level; hence, individual curves can be fit with common fixed 
effects. In the figure, each column represents the case of one 
subject and each row represents the case of one area. For 
example, column 1 from the left presents subject l's source 
strength-time curves in the 5 areas with the subjectl/areal 
being the curve for area 1 (left Ml) and subject l/area5 for 
area 5 (PM&SMA). It was noted that, although the 8 subjects 
had different curve shapes, the general patterns of the curves 
were similar and the peaks of the curves located consistently 
near the force onset time (-time 0) for all the cases. 

DISCUSSION 

This study found that source strength (measured by 
current density [CD]) of scalp EEG signals during voluntary 
handgrip was time-dependent. The source increased almost 
linearly from about 1000ms before and peaked at the force 
onset in a number of motor function-related cortical areas; 
the patterns of the time-dependent variation of the CD were 
similar among the cortical regions. 
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Table 1. Analysis of Covariance Based on the Mixed-Effects Polynomial Regression Model. The "Time" Variable was Normalized 
into [0,1] in the Analysis. The p-Values Based on ANOVA F Tests are Presented for Global Fixed Effects Over all 
Subjects 
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Fig. (3). The original data and fitted curves for each area within each subject from the mixed-effects polynomial regression model. Each 
column represents the case of one subject and each row represents the case of one area. Area 1, area2, area3, area 4, and area5 denoted the 
left Ml, right Ml, left SI, right SI, and PM and SMA, respectively. For example, the sub-plot subject5/area3 denotes the case of subject 5 
with area3: Left SI. 



Time-Dependent Source Strength Variation in Control- 
ling a Voluntary Handgrip 

Our results indicated that from as early as 1 s before the 
force onset the source strength had elevated in cortical 
regions examined (Fig. 1, 2 and 3) and peaked near the time 
of force onset. Cortical activation at early times may be more 
related to attention and general preparation of a voluntary 
motor action [23]. At early stages of motor preparation, 



perhaps only a limited number of neurons are involved in the 

task, and the activation level might not critically depend on 
the force exerted [24-29]. As the time approaches motor task 
initiation, more activities related to specific programming 
and execution of the task occur and this may require prog- 
ressively more cortical neurons to participate and increase 
their activation levels. This postulation is supported by the 
almost linear increases of the source strength from -1000ms 
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to the force onset (time 0, Fig. 2). The observation of the 
source strength at the highest level around the force onset 
could be contributed by participation of neurons that 
executed the handgrip and neurons whose role involved in 
modulating the grip force. Cheney [30] suggested that about 
100ms before the EMG onset, force-magnitude modulation, 
known as increasing discharge rate and/or recruitment of 
more neurons, might have been established. About 100ms or 
more after the EMG onset, this modulation is further 
adjusted so that the cortical activity increases according to 
the level of force generated by the muscle. In summary, the 
high levels of cortical activation during the time window of 
± 100ms movement onset are highly likely to be related to 
specific motor task programming, execution, force modula- 
tion, and perhaps as well as multi-level sensory feedback 
reaching to the motor cortices [30]. 

During the later phase of the holding period (~200ms 
after the force onset), the source strength began to decline in 
all the analyzed areas (Fig. 2, Fig. 3). The reduction in the 
cortical activity during the holding phase could be a result of 
ceasing activation of the neurons involved in motor action 
programming and execution. Many cortical cells increase 
discharge rate before the force onset but most of these 
neurons were silent during the period of holding the force. 
For example, more than 60% of the recorded cells in monkey 
SMA and cingulated cortex were phasically active before 
onset of finger force but only about 10% of the neurons were 
tonically active during the holding phase of the force task 
[31]. In addition, motor neurons in the spinal cord adapt to 
lower firing rate after initial high-frequency activation [32] 
and this adaptation could also occur in cortical neurons [33], 
which may as well explain weakening of the source strength 
during later phase of the force holding. 

Similar Source Strength Variation Time Courses among 
the Cortical Regions 

All the examined regions showed time-dependent source 
strength changes during the control process of the muscle 
contraction. The left Ml and SI exhibited stronger current 
density (CD) during almost the entire time course (Fig. 2), 
which was expected because the handgrip occurred on the 
right side. An interesting observation was that although the 
level of the CD seemed to differ among the cortical areas, 
the pattern of the source strength variation in the time course 
was very similar, that is all the areas exhibited a reversed 
"V-shape" CD time course (Fig. 2). This was confirmed by 
the statistical F-test that all terms of polynomial coefficients 
for the nonlinear current density time course were significant 
while the ANOVA test for the factor of "area" was not. 
Previous animal studies involving single-cell recording 
reported that the Ml is primarily for controlling movement 
execution variables, such as the amount of force to be 
exerted by the muscle [1, 30] and direction of the force 
production [34]. The Ml also encodes higher functions, such 
as movement programming and preparation [35]. The SI is 
known to participate in handgrip force control [36] and 
provide sensory information to Ml during voluntary motor 
actions [37]. The SMA has been shown to be important in 
programming sequential finger movements [38], storing 
information necessary for the orderly performance of multi- 
ple movements of the arm [39], and bimanual coordination 
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[40, 41]. The PM is considered to contribute to the selection 
of motor actions on the basis of visual cues [42, 43]. The 

monotonic relationship between cell activity and muscle 
output holds for different areas in the motor cortex in a series 
of single-cell animal studies with slightly different experi- 
ment paradigms [1, 44-46]. Specifically, the Ml and PM 
cortices share the same force-firing rate pattern [46]. Two 
populations of neurons (phasic-tonic firing and tonic 
decreasing firing) in the Ml and PM cortex have been found 
to be related, either positively or negatively, to force output 
in a more recent study [47]. The relationship between the 
pyramidal tract neurons and force output was monotonic 
over only a small portion of the force range [45]. In 
summary, all the examined areas play important motor 
function, which is underscored by their contributions to the 
cortricospinal tract (CST), a fundamental pathway for 
voluntary actions of extremity muscles (one third of the CST 
is contributed by the Ml , one third by the PM and SMA, and 
the remaining one third by the SI [48]). 

One limitation of the single-cell recording in animal 
models, however, is the difficulty to investigate activities of 
multiple cells in multiple cortical regions concurrently. Thus, 
it is not well known, from the animal data, the integrated 
activation patterns of multiple brain areas. This limitation 
can be overcome by modern neuroimaging techniques and 
abundant observations have suggested high degrees simi- 
larities in activation patterns of many sensorimotor function- 
related cortical areas. For example, neuroimaging data have 
shown a proportional relationship between the levels of 
muscle and activation, and this relation is similar among 
many cortical regions [6, 36, 49, 50]. A number of studies 
have reported muscle fatigue-related brain activation 
alterations and a prominent finding by these studies was that 
the cortical signal changing pattern was very much alike 
among many sensorimotor regions [4, 5, 17]. The human 
neuroimaging data support the observation of similar source 
strength time courses among the cortical areas in the current 
study and suggest synchronized parallel functions among 
multiple sensorimotor centers. 

Compared with previous findings of EEG source signals 
generated during higher levels of handgrip contractions, the 
relative source strength (CD) of the SMA and PM (relative 
to CD of Ml and SI) was greater for the 20% MVC of the 
current study than relative source strength of the SMA and 
PM at higher levels of handgrip [51]. It is not clear why this 
was so but it could be that fewer cells in the primary 
sensorimotor areas need to be activated to plan, execute and 
maintain a low level of muscle contraction but the level of 
activation in the SMA and PM may not change as much 
across degrees of muscle contraction. 

In conclusion, the results show that the level of cortical 
activity in controlling a voluntary motor action is time- 
dependent with the highest activity occurs around the action 
onset. The low-to-high-to-low source strength pattern from 
the preparation to execution and to holding phases is similar 
among several major sensorimotor areas, suggesting syn- 
chronized parallel functions among them. The findings also 
indicate that the current density reconstruction technique 
based on surface EEG recorded during voluntary motor 
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activities has adequate time and spatial resolution to detect 
the time-dependent source strength in various brain regions. 



ACKNOWLEDGEMENTS 

This work was supported by research funds of 
Department of Physical Medicine and Rehabilitation at the 
Cleveland Clinic 

REFERENCES 

[I] Evarts EV. Relation of pyramidal tract activity to force exerted 
during voluntary movement. J Neurophysiol 1968; 3 1 : 14-27. 

[2] Conrad B, Wiesendanger M, Matsunami K, Brooks VB. Precentral 
unit activity related to control of arm movements. Exp Brain Res 
1977; 29: 85-95. 

[3] Fetz EE, Cheney PD. Postspike facilitation of forelimb muscle 
activity by primate corticomotoneuronal cells. J Neurophysiol 
1980;44:751-72. 

[4] Liu JZ, Shan ZY, Zhang LD, Sahgal V, Brown RW, Yue GH. 

Human brain activation during sustained and intermittent 

submaximal fatigue muscle contractions: an fMRI study. J 

Neurophysiol 2003; 90: 300-12. 
[5] van Duinen H, Renken R, Maurits N, Zijdewind I. Effects of motor 

fatigue on human brain activity, an fMRI study. Neuroimage 2007; 

35: 1438-49. 

[6] Dettmers C, Fink GR, Lemon RN, et al . Relation between cerebral 
activity and force in the motor areas of the human brain. J 
Neurophysiol 1995; 74: 802-15. 

[7] Siemionow V, Yue GH, Ranganathan VK, Liu JZ, Sahgal V. 
Relationship between motor activity-related cortical potential and 
voluntary muscle activation. Exp Brain Res 2000; 133: 303-1 1. 

[8] Tecchio F, Zappasodi F, Melgari JM, Porcaro C, Cassetta E, 
Rossini PM. Sensory-motor interaction in primary hand cortical 
areas: a magnetoencephalography assessment. Neuroscience 2006; 
141(1): 533-42. 

[9] Michel CM, Murray MM, Lantz G, Gonzalez S, Spinelli L, Grave 
de Peralta R. EEG source imaging. Clin Neurophysiol 2004; 115: 
2195-222. 

[10] Babiloni F, Babiloni C, Locche L, Cincotti F, Rossini PM, 
Carducci F. High-resolution electro-encephalogram: source 
estimates of Laplacian-transformed somatosensory-evoked 
potentials using a realistic subject head model constructed from 
magnetic resonance images. Med Biol Eng Comput 2000; 38: 512- 
19. 

[II] Yao J, Dewald JP. Evaluation of different cortical source 
localization methods using simulated and experimental EEG data. 
Neuroimage 2005; 25: 369-82. 

[12] Winterer G, Mulert C, Mientus S, et al. P300 and LORETA: 

comparison of normal subjects and schizophrenic patients. Brain 

Topogr 2001; 13:299-13. 
[13] Pae JS, Kwon JS, Youn T, et al. LORETA imaging of P300 in 

schizophrenia with individual MRI and 128-channel EEG. 

Neuroimage 2003; 20: 1552-60. 
[14] Lantz G, Michel CM, Pascual-Marqui RD, et al. Extracranial 

localization of intracranial interictal epileptiform activity using 

LORETA (low resolution electromagnetic tomography). 

Electroencephalogr Clin Neurophysiol 1997; 102: 414-22. 
[15] Worrell GA, Lagerlund TD, Sharbrough FW, et al. Localization of 

the epileptic focus by low-resolution electromagnetic tomography 

in patients with a lesion demonstrated by MRI. Brain Topogr 2000; 

12:273-82. 

[16] Thut G, Hauert CA, Morand S, Seeck M, Landis T, Michel C. 
Evidence for interhemispheric motor-level transfer in a simple 
reaction time task: an EEG study. Exp Brain Res 1999; 128: 256- 
61. 

[17] Liu JZ, Dai TH, Sahgal V, Brown RW, Yue GH. Non-linear 
cortical modulation of muscle fatigue: A functional MRI study. 
Brain Research 2002; 957: 320-29. 

[18] Fuchs M, Drenckhahn R, Wischmann HA, Wagner M. An 
improved boundary element method for realistic volume-conductor 
modeling. IEEE Trans Biomed Eng 1998; 45: 980-97. 



[19] Fuchs M, Kastner J, Wagner M, Hawes S, Ebersole JS. A 

standardized boundary element method volume conductor model. 

Clin Neurophysiol 2002; 113: 702-12. 
[20] Hedeker D, Gibbons RD. Longitudinal data analysis. Chapter 5. 

John Wiley & Sons 2006. 
[21] Davis CS. Statistical methods for the analysis of repeated 

measurements. New York: Springer 2002. 
[22] Verbeke G, Molenberghs G. Linear mixed models for longitudinal 

data. New York: Springer 2000. 
[23] Ball T, Schreiber A, Feige B, Wagner M, Lucking CH, Kristeva- 

Feige R. The role of higher-order motor areas in voluntary 

movement as revealed by high-resolution EEG and fMRI. 

Neuroimage 1999; 10: 682-94. 
[24] Buys EJ, Lemon RN, Mantel GW, Muir RB. Selective facilitation 

of different hand muscles by single corticospinal neurones in the 

conscious monkey. J Physiol 1986; 381 : 529-49. 
[25] Ashe J. Force and the motor cortex. Behav Brain Res 1997; 87: 

255-69. 

[26] Shen L, Alexander GE. Neural correlates of a spatial sensory-to- 
motor transformation in primary motor cortex. J Neurophysiol 
1997;77:1171-94. 

[27] Sakata H, Taira M, Kusunoki M, Murata A, Tanaka Y. The TINS 
Lecture. The parietal association cortex in depth perception and 
visual control of hand action. Trends Neurosci 1997; 20: 350-7. 

[28] Murata A, Gallese V, Luppino G, Kaseda M, Sakata H. Selectivity 
for the shape, size, and orientation of objects for grasping in 
neurons of monkey parietal area AIP. J Neurophysiol 2000; 83: 
2580-601. 

[29] Morrow MM, Miller LE. Prediction of muscle activity by 

populations of sequentially recorded primary motor cortex neurons. 

J Neurophysiol 2003; 89: 2279-88. 
[30] Cheney PD. Role of cerebral cortex in voluntary movements. A 

review Phys Ther 1985; 65: 624-35. 
[31] Cadoret G, Smith AM. Comparison of the neuronal activity in the 

SMA and the ventral cingulate cortex during prehension in the 

monkey. J Neurophysiol 1997; 77: 153-66 
[32] Kernell D, Monster AW. Motoneurone properties and motor 

fatigue: an intracellular study of gastrocnemius motoneurones of 

the cat. Exp Brain Res 1982; 46: 197-204. 
[33] Hepp-Reymond M-C, Wannier TMJ, Maier MA, Rufener EA. 

Sensorimotor cortical control of isometric force in the monkey. 

Progress in Brain Res 1989; 80: 451-63. 
[34] Georgopoulos AP, Schwartz AB, Kettner RE. Neuronal population 

coding of movement direction. Science 1986; 233: 1416-9. 
[35] Michelon P, Vettel JM, Zacks JM. Lateral somatotopic 

organization during imagined and prepared movements. J 

Neurophysiol 2006; 95: 81 1-22. 
[36] Dai TH, Liu JZ, Sahgal V, Brown RW, Yue GH. Relationship 

between muscle output and functional MRI-measured brain 

activation. Exp Brain Res 2001; 140: 290-300. 
[37] Berne RM, Levy MN, Koeppen BM, Stanton BA. Physiology-5 th 

edition. Mosby 

[38] Roland PE, Larsen B, Lassen NA, Skinhoj E. Supplementary motor 
area and other cortical areas in organization of voluntary 
movements in man. J Neurophysiol 1980; 43: 1 18-36. 

[39] Tanji J, Shima K. Role for supplementary motor area cells in 
planning several movements ahead. Nature 1994; 371: 413-6. 

[40] Brinkman C. Lesions in supplementary motor area interfere with a 
monkey's performance of a bimanual coordination task. Neurosci 
Lett 2003; 27: 267-70. 

[41] Brinkman C. Supplementary motor area of the monkey's cerebral 
cortex: short- and long-term deficits after unilateral ablation and 
the effects of subsequent callosal section. J Neurosci 1984; 4: 918- 
29. 

[42] Halsband U, Passingham RE. Premotor cortex and the conditions 
for movement in monkeys (Macaca fascicularis). Behav Brain Res 
1985;18:269-77. 

[43] Passingham RE. Cues for movement in monkeys (Macaca mulatta) 
with lesions in premotor cortex. Behav Neurosci 1986; 100: 695- 
703. 

[44] Hepp-Reymond MC, Wyss UR and Anner R. Neuronal coding of 
static force in the primate motor cortex. J Physiol (Paris) 1978; 74: 
287-91. 

[45] Cheney PD, Fetz EE. Functional classes of primate 
corticomotoneuronal cells and their relation to active force. J 
Neurophysiol 1980; 44: 773-91. 



Cortical Activation in Handgrip Contraction 

[46] Hepp-Reymond M, Kirkpatrick-Tanner M, Cabernet L, Qi HX, 
Weber B. Context-dependent force coding in motor and premotor 
cortical areas. Exp Brain Res 1999; 128: 123-33. 

[47] Hepp-Reymond MC, Hiisler EJ, Maier MA, Ql HX. Force-related 
neuronal activity in two regions of the primate ventral premotor 
cortex. Can J Physiol Pharmacol 1994; 72: 571-9. 

[48] Ghez C, Krakauer J. The Organization of Movement. In: E.R. 
Kandel, J.H. Schwartz and T.M. Jessell (Eds.), Principles of Neural 
Science (4th ed). McGraw-Hill, New York, 2000; 653-73. 



The Open Neuroimaging Journal, 2011, Volume 5 239 

[49] Thickbroom GW, Phillips BA, Morris I, Byrnes ML, Mastaglia FL. 

Isometric force-related activity in sensorimotor cortex measured 

with functional MRI. Exp Brain Res 1998; 121: 59-64. 
[50] Cramer SC, Weisskoff RM, Schaechter JD, et al. Motor cortex 

activation is related to force of squeezing. Hum Brain Mapp 2002; 

16: 197-205. 

[51] Wang XF, Yang Q, Fan Z, Sun CK, Yue GH. Assessing time- 
dependent association between scalp EEG and muscle activation: A 
functional random-effects model approach. J Neurosci Methods 
2009; 177: 232-40. 



Received: November 06, 2010 Revised: May 10, 201 1 Accepted: June 27, 201 1 

© Yang et al.; Licensee Bentham Open. 

This is an open access article licensed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by- 
nc/3.0/), which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited. 



